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PREFACE 


This  report  describes  work  done  by  Or  B.  0.  Sivazlian  and  Dr  J.  F. 
Mahoney,  Department  of  Industrial  and  Systems  Engineering,  the  University  of 
Florida,  Gainesville,  Florida  32611  under  Contract  No.  F08635-83-C-02D2 
sponsored  by  the  Air  Force  Armament  Laboratory  (AFATl).  The  program  manager 
was  Mr  Daniel  A.  Mclnnis  (DLYW). 

The  work  was  initiated  under  a  1982  USAF-SCEFF.  Summer  Faculty  Research 
Program  sponsored  by  the  Air  Force  Office  of  Scientific  Research  conducted  by 
the  Southeastern  Center  for  Electrical  Engineering  Education  (SCEFE)  under 
Contract  No.  F49620-82-C-0035. 

This  work  addresses  itself  to  the  problem  of  computing  the  uncertainty 

associated  with  the  probability  of  kill  due  to  blast  in  the  presence  of 

aiming  error  and  in  the  absence  of  fragmentation.  Let  Pk  be  the  probability 

of  kill  due  to  blast  alone.  The  assumption  is  made  that  (1)  Pk  is  unity 

between  the  center  of  the  blast  and  a  distance  A  from  the  center,  (2)  Pk  is 

negligible  beyond  a  certain  distance  B  from  the  center  (B  >  A),  (3)  Pk 

decreases  linearly  between  A  and  B,  (4)  the  aiming  errors  in  the  direction  of 

range  and  deflection  are  unbiased,  independent  of  each  other  having  a  Gaussian 

distribution  with  respective  standard  deviations  a  and  o  .  It  is  shown  that 

x  y 

Pkp  can  be  expressed  mathematically  by  integrals  involving  the  modified  Bessel 
function  of  the  first  kind  I q ( - ) .  Moreover,  under  the  assumption  that  the 
input  parameters  A,  B,  a  and  a  are  not  known  with  certainty,  but  are 

*  y 

■  haracteri zed  by  their  first  two  moments,  explicit  approximate  expressions  are 
obtained  for  E F. P k and  Var[Pkh], 


i 


The  author  has  benefited  from  helpful  discussions  with  several  people. 
Particular  thanks  are  due  to  Mr  Jerry  P.  Bass,  Mr  Daniel  A.  Mclnnis  and  Mr 
Charles  A.  Reynolds,  all  from  DLYW,  who  have  read  the  report  and  have 
contributed  to  it  through  helpful  comments. 

The  report  is  the  fourth  of  a  series  dealing  with  the  uncertainty 
associated  with  various  weapon  effectiveness  indices  and  details  methodologies 
and  techniques  used  in  computing  such  uncertainties  in  the  presence  of  error 
in  the  ’nput  parameters. 

The  Public  Affairs  Office  has  reviewed  this  report,  and  it  is  releasable 
to  the  National  Technical  Information  Service  (NTIS),  where  it  will  be 
available  to  the  general  public,  including  foreign  nationals. 
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SECTION  I 


INTRODUCTION 

This  report  considers  the  problem  of  estimating  the  probability  of  kill 
due  to  blast,  Pkb,  in  the  presence  of  aiming  error.  Prior  to  solving  the 
estimation  problem,  it  is  necessary  to  compute  a  mathematical  expression  for 
Pkb  in  order  to  apply  the  usual  statistical  techniques  to  arrive  at  confidence 
intervals  for  Pkb.  In  Section  II,  the  derivation  of  Pkb  is  carried  out  Rased 
on  several  explicitly  stated  assumptions.  In  Section  III,  the  estimation 
problem  for  the  case  when  range  aiming  error  equals  deflection  aiming  er-or  is 
solved.  In  Section  IV,  the  general  estimation  problem  is  accomplished  w^en 
these  two  aiming  errors  are  unequal.  Finally,  Section  V  provides  some 
conclusive  remarks. 

No  particular  problem  is  encountered  in  Section  III  when  the  two  aiming 
errors  are  equal.  It  will  be  seen  that  Pkb  is  expressible  in  terms  of  the 
error  function.  The  estimates  E[Pkbl  and  Var[Pkb]  can  be  obtained  in  closed 
form,  and  again  will  involve  the  error  function.  However,  it  will  be  seen 
that  the  most  general  expression  for  Pkb»  when  {  a^,  involves  integrals  of 
the  modified  Ressel  function  of  zero  order  Ib(*l.  Although  this  result  's 
interesting  in  itself,  nevertheless,  it  becomes  difficult  to  overcome  * "o 
mathematical  intricacies  introduced  by  the  presence  of  this  function.  To 
arrive  at  closed  form  expressions  for  the  mean  and  variance  of  nkb,  in  the 
general  case,  it  will  be  necessary  to  use  approximation  techniques  which 
result  in  rather  lengthy  expressions  for  the  estimators. 


SECTION  II 


MATHEMATICAL  MODEL  FOR  Pkb 


1.  Background 

In  this  section,  a  mathematical  model  for  the  probability  of  kill  due  to 
blast,  Pkb,  in  the  presence  of  aiming  error  is  developed.  The  basic  situation 
that  one  is  facing  consists  of  the  following. 

A  wpapon  whose  main  effect  is  kill  due  to  blast  is  delivered  from  air  to 
a  targe*1,  point  located  on  the  ground  surface.  The  target's  position  is 
stationary.  The  weapon  may  not  hit  directly  the  target  due  to  the  presence  of 
aiming  errors.  These  errors  are  assumed  to  be  unbiased;  that  is,  centered  at 
the  location  of  the  target.  The  target  may  or  may  not  bp  killed  by  the  effect 
of  blast  depending  on  whether  the  weapon  impacts  close  to  or  away  from  the 
target.  The  probability  of  kill  due  to  blast  is  related  to  the  distance 
between  weapon  and  target  by  a  well-defined  mathematical  function.  In 
addition,  the  aiming  error  is  not  known  precisely  but  is  expressed  by  a 
probability  density  function  which  provides  a  mathematical  formula  for 
computing  the  probability  that  the  weapon  will  impact  in  an  interval  du  dv 
close  fe  a  point.  (u,vl  on  the  ground  surface. 

The  technique  that  will  be  used  to  compute  the  expression  for  the 
probability  of  kill  due  to  blast,  Pkb,  in  the  presence  of  aiming  error  is 
based  on  the  laws  of  conditional  probability.  Ultimate!  •v-  Pkb  is  not  going  to 
depend  >n  the  position  of  the  target  if  the  aiming  error  is  unbiased,  and  if 
weapon  delivery  can  theoretically  result  in  a  point  of  impact  which  can  be 
anywher.'  .on  the  ground  surface.  On  the  other  hand,  nkb  will  depend  on: 

a.  1  he  parameters  specifying  the  functional  form  relatiny  probability  of 
kill  to  distance . 

b.  The  statistical  parameters  of  the  aiming  error  distribution. 


2. 


Assumptions 

The  following  assumptions  pertaining  to  this  situation  will  be  made: 

a.  Both  the  target  and  the  weapon  are  idealized  as  points,  and  the 
weapon  is  aimed  at  the  target. 

b.  The  direction  of  the  weapon  delivery  range  and  deflection  are, 
respectively,  parallel  to  the  x  and  y  coordinate  axes  on  the  ground  plane. 
Since  the  coordinate  system  can  be  arbitrarily  selected,  there  is  no  loss  in 
generality  in  making  this  specific  assumption.  The  position  of  the  target  has 
coord i nates  ( x  ,y )  . 

c.  The  aiming  error  (distance)  in  each  of  the  x  and  y  directions  are 

independently  and  normally  distributed  with  respective  means  x  and  y  and 

standard  deviation  o  and  o  .  Let  (du  dv)  be  the  infinitesimal  rectangle, 
x  y 

close  to  the  point  (u,v)  at  which  the  weapon  impacts,  and  define  the  random 
variables  U  and  V  which  measure,  respectively,  the  distances  between  the 
target  point  and  the  weapon  impact  point  along  the  abscissa  and  the 
ordinate.  Then,  the  probability  that  the  weapon  will  impact  in  the  rectangle 
(du  dv)  is: 

f,,  v  ( u-x ,  v-y)  du  dv  =  „  *  ;  exp[  -  * 

’  '  '  x  y  2o  x 

exp[-  ^v~j^  ■]  du  dv  .  (1) 

2<j 

y 

d.  The  probability  of  kill  due  to  blast  at  a  point  (x,y)  given  that  the 
weapon  impacts  at  (u,v)  depends  only  on  the  distance  r  between  (x,y)  and  <  \i  ,v) 

on 

r  =  /  (x-u)'’*  +  (y-v)? 


and  is  given  by 


1 


0  <  r  <  A 


Vr) 


8-r 

B-A 

0 


A  <  r  <  R  . 
B  <  r  <  » 


(?) 


Thus,  if  the  target  is  located  at  a  distance  less  than  A  from  the  center  of 
the  blast,  the  probability  of  kill  is  unity.  If  the  target  is  located  at  a 
distance  exceeding  B,  (R  >  A),  the  damage  to  the  target  is  negligible  and  the 
probability  of  kill  is  zero.  Between  A  and  3,  the  probability  of  kill  is  a 
linear  function  of  the  distance,  if  the  delivery  of  the  weapon  results  in  a 
direct  hit  on  the  target  (r=0),  tne  probability  of  kill  is  unity. 

e.  Fragmentation  effect  is  neglected.  This  implies  that  the  target  is 
not  fragment  sensitive.  Or,  if  weapon  fragmentation  exists,  its  effect  on  the 
target  is  negligible,  hence,  not  resulting  in  a  kill. 

f.  In  general,  the  point  in  space  from  which  weapon  is  delivered  is 
nonstationary,  the  weapon  is  subject  to  ballistic  errors,  etc.  It  shall  be 
assume!  that  all  these  factors  combine  into  a  single  source  of  error  which  is 
i "co-p ’rated  in  the  aiming  error. 

g.  The  blast  does  not  contribute  to  the  aiming  error. 

3.  Expression  for  > 


The  probability  of  kill  due  to  blast  is 


oo  » 

P..h  =  /  f  [Probability  of  kill  at  (x,y)  ]  weapon  impacts  at  (u,v)l  X 

„oo  _uO 


[Probability  that  the  weapon  impacts  between  (u,v)  and 
(u+du,  v+dv)] 


=  I 


UO  ■—  — —  —  —  —  —  —  —  -  —  —  ■- 

f  p,[/  (x-u)7  +  (y-v)^] 


f  u(u-x,  v-y)  du  dv 


=  /  /  p.  [/ (x-u)2  +  (y-v) 2 ] - - -  Lmu. 


_aO  —  00 


2"Vye<Pl_  *2 


exp[  ■ 


( v-v)' 


2  o 


2 


du  dv  .  (3) 


In  the  (x,y)  plane,  translate  the  axis  in  such  a  way  that  the  origin  coincides 
with  (x,y),  which  is  the  location  of  the  target.  Then  letting 


X  =  u-x  and  Y  =  v-y 


one  obtains 


pkb  =  /  /  Pk  t/  X2  +  Y2)  27^  exp[-  (-4  +  -4)]  dx  dY  • 


x  y 


2j  2o 

x  y 


Changing  now  to  polar  coordinates,  let  X=rcos  6  and  Y=rsin  e.  Inis  yields 


2tt 


2  2„  2.2, 


kb 


-I  /  Pk(r)?^exp[.(^^+^LnO]]  rdrdJ 


0  n  kV’7  2"  Vy 


?o 


2  o 


(5) 


Substituting  for  the  value  of  Pk(r)  as  given  in  (2),  one  obtains 


kb  2w  a  a 

x  y  0  0 


?tt  A  2  2  2.2 

r  cos  -j  r  sin  •  ■ 

J  ,  expt-L - + - j—Ji  r  nr  -is 

2a^  2o4 

x  y 


2tt  B 


j.  r  c  r  B-f\  r  ci  i_us  o  r  sin  tin 

irrrf,  /.  (rod  expi-(— rr-*— p— ii  r 


2  2a  2  .  2. 
r  cos  0  .  r  sin  9- 


dr  de 


x  y  0  A 


2o 


2a . 


(6) 


Notice  here  that  one  is  dealing  with  integrals  of  the  form 


2n  2  ?  2  .  9 

,  r  r  cos  a  r  sn  a ,  -  . , 

K  r.  =  j  exp  - ^ —  a  — — 5 — Jj  do  . 


0 


2a 


2a 


(/) 


In  Appendix  A  it  is  shown  that 


where 


and  !  • 

hecorr>s 


K(r )  =  p'3r  In(hr2^ 


(«' 


.  -  1  il  + 

rl  ~  ~r  i  — tt-  * 


l  , 

T  T  TJ 

a  a 

x  y 


(9) 


.  _  1  ,1  1  , 
h  ‘  T  'T  "  “T 

o  o 

y  x 


(10) 


1  is  a  modified  Ressel  function  of  zero  order.  Expression  (6)  for  P, 


.  A  2  „ 

P,.^  =  -7-7-  j  r  e"ar  I0(br  )  dr 


k .  t>  oo  n 

x  y  0 


+  J  r  P  ^  I0(br?)  dr  • 

x  v  A 


an 
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section  III 


ESTIMATION  OF  E[?kb]  AND  Var[Pkb]  WHEN  =  a 

1 .  Background 

In  this  section,  the  estimation  problem  for  Pkb  will  be  investigated  for 
the  particular  case  when  the  range  aiming  error  is  equal  to  the  deflection 
aiming  error  (a^  =  °y  =  °)*  This  case  is  of  practical  significance  in 
guided  weapons  for  which  the  aiming  error  in  all  directions  is  about  the 
same.  Moreover,  the  problem  itself  allows  one  to  obtain  closed  form 
expressions  for  the  estimates  E[Pkb]  and  Var[Pkb]  involving  the  error  functioi 

2 

w 

1  ?  -  o 

$(z)  =  -  f  e  dw.  (12) 

✓~F  n 


Now,  since  lg(0)  =  1,  it  follows  from  (11)  that 

2 


2 


,  1  /  2a  ^  1  fB  ,-B-r-,  2 o?  , 

’kb  3  T  /  r  R  dr  +  “2  !  ^  r  e  f1r  ' 

a  0  o  A 


This  expression  can  be  reduced  into  three  integrals  as  follows: 


,  A  ‘  2  „  .  B  '  ,  2 

Pkb  »  ^2  /  e  rdr.X.Lj  r  dr 

a  0  o  A 


1  1  2  2?  , 
iw\  *  T  1  r  e  dr 

a  A 


A  ~  }?  M  B 

*  /  -  (1IP  )  +  ITT  J  'dl< 

0  ‘  A 


2a 


)  +  -R--A-  J  r  <V 


Noting  that  the  two  first  integrands  are  exact  differentials  and  integrating 
by  parts  the  last  integral,  one  obtains 


which  upon  simplification  results  in 


B  "  — 7 

pkl>  -  1  -  iPT  /  e  2°  dr- 


Let  x  -  r/o,  then 


B  2 
—  x 
a  -  T~ 

=  I  -  —  /  e  ‘  dx 


kb 


R-A 


A 

o 


which  mav  be  written  as 


kb 


1  - 


v  a  r  1 


B 
a 


(B-A) 


[_I_  I  e  2  dx  -  4=  /  e  2  dx] 


/  2j  n 


/  2ir  0 


Using  (IP)  one  finally  obtains 


kb 


1  /  2n  J  r  .  /•  B  \ 

1  -  — —  L*  (-j 


(R-A) 


•  $]■ 


(13) 


(14) 


It  will  In*  implicitly  assumed  that  the  parameters  A,  R,  and  a  are  not  known 
exactly,  hut  are  subjected  to  estimation  error.  Further,  A,  B,  and  o  will  be 


v-1;  nr.  v_  r  irj  L  . 


assumed  to  be  mutually  independent  and  defined  by  their  respective  means 
A,  B,  a  and  their  respective  variances  Var[A],  Var[B],  and  Var[o]. 

Recalling  that  Pkb  in  (14)  is  a  function  of  the  three  parameters  A,  R, 
and  a,  one  may  write  (14)  as 


pkb(M,a)  =  1 


T^f  I‘(f)  -•(£)] 


(15) 


Let  A,  B  and  o  refer,  respectively,  to  the  mean  of  A,  B  and  o.  Expanding  Pkb 
about  the  point  (A,B,a)  one  obtains  up  to  the  first  order  terms 


-  3P 


Pkb(A,Bto)  =  Pkh(A,B,o)  +  (A-A)  ^ 

-  3P 


kb 


A,B,a 


♦  (B-B) 


kb 


3B 


A,R,o 


3P,, 

/  —  kb 

+  (o-o)  -37- 


'16' 


A  ,R  ,0 


2.  Estimation  of  E[Pk(j] 

Taking  expectations  on  both  sides  of  (16)  yields  as  a  first  approximation 


Erpkb(A,R,o)l  -  pkb(A,B,a) 


'17) 


or  using  (15) 


■  1  - 


2ti  0  r  *  r~K' 1 


__  [♦  0  -  ♦  0] 

(B-A)  a  a 


Let 


Then 


u( a,b,u)  =  inf)  -  H^)]  • 


(1«) 


(IP) 


E[Pkbl  *  1  -  0  U  (A.B.a)  . 


(20) 
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3. 


Estimation  of  VarCP^] 


First,  expression  (16)  is  written  as 


PkblA,B,a)  -  Pkb(A,B,o)  =  (A-A) 


apkb 

3A 

A,B 

-  3P.  , 

+  (B-B)  -$■ 


-  3P 


-  +  (o-o) 


kb 


A,B  ,o 


3  a 


A  ,R  ,a 


(21 


It  was  assumed  that  each  of  the  random  variables  A,  B,  and  o  is  mutually 
i ndependent .  Thus,  squaring  and  taking  the  variance  on  both  sides  of  (21) 
results  in: 


Varfl1 


kb' 


3 P,  .  ?  3P  ? 

(-#) _ Var[A]  +  K-f ) _ Varpl] 


3P.. 

+  (-A 

^  3o  ; 


A,B,o 
2 


A,B  ,o 


"A.B.o 


Var[a]  . 


(22 


ro  o  fa-'-  the  partial  derivatives,  form  (13)  of  °kb  will  be  used  rather  than 


form  (1J).  First  compute  3Pkby^A: 


.2 


3P 


kb 


1 


( B-A)  A 


B  _  2  ,  o 

r  ?o  .  1  2o 

f  e  ’  drt^e 


2 


o  /  2ti  r  1 


B  2 
—  x 
0 


A  2 
—  x 
a 


e  ?dx  -  J  p  'cix] 


1  ?.< 

+  — r  e 


( B-A)  /2«  n 


✓  2  n  n 
,2 


B-A 


J  '  ?  t*  (J)  -  *£)]  ^  ^  *  ?0<' 


(B-A) 

Using  the  definition  (Id)  yields  finally 


l 


3P ..  ,  -  -2 

~Tir  ~  ~  aZK  l'( A,B ,0)  +  r_a  e 


(23 


Next  3P|<^3R  is  computed: 


3P 


kb 


1 


38 


(B-A )  A 


f  20  w 

j  e  dr 


XT 


_R 

2o 


R  2  A  2 

^  i_!_ /V  **  *  . 

(R-Ap  ST*  0  /7S"  n 


Using  (19),  this  expression  nay  he  written  as 


3P 


kb  a  /7n  -P.  a,  A',|  1 

=  - 7  L*l-J  -  n-j]  -  5T 


R2 

77 


3B  (B-A 7  ‘ 


R-A 


Finally,  dP^^/do  is  computed.  Using  (13)  yields 

2 


3P,  .  ,82  T  0  2 

kb  _  1  r  r  ,  ,,  -3  2 a 

To - -7("2)°  e 

A  2 

r 


dr 


3P 


kb_ 

a 


1 

(B -A )"'  a 


1 

JWTJT 


R  “  T7 

/  -r  (-£•)  e  °  dr 
A  o  ? 

/8  r  d  (e  2?)  . 


Integrating  by  part  results  in  (see  Section  1 1 1  - 1 ) : 


3P 


kb 


3cr 


1 

(B'-A)'a 


"  rr  _  7T  b  '  7T 

[Re  a-Ae  a-/e  °drj 


,2 


B 


1  ,  «  20"  D  77; 

-  [A  e  “ B  e  J 

B  2  A  2 

_  -  _  x  -  x 

/2  it  ;  1  2  .  1  ”7  .  1 

+  J  e  dx  -  —  J  e  dxj. 


m  o 


/7n  0 


1 

XT  p 


Using  the  definition  (19)  gives 


11 


Ml 


(r-aT  o  e 


_  f.u  a  T  C  '  r  .  'M  r\  a  ' 

-Be  1  +  -rTa  *  ^o;l  ‘ 


(25) 


Substituting  (23),  (24),  and  (25)  in  (22)  yields 

—  “  -2 

Var[P,.  1  =  [-  -2-  ll(A,B,o)  +  e  ?0  ]?  Var[Al 
KD  B-A  R-A 

’  ~=? 

+  [~  tl( A ,B ,o )  -  -i-  e  7a  }2  VarfB] 
B-A  B-A 


"  ~  —  0 

+  [ - - -  ( A  e  ?a  -  B  e  2°  )  +  ')(A,B,ol]?  Varr^l.  (?fi' 

(B-7f)  a 


4.  Example 

A  guided  weapon  whoso  main  effect  is  blast  is  aimed  at.  the  renter  n  of  a 

target.  Because  of  error,  the  weapon  will  not  necessarily  impact  at  0.  This 

error,  measured  as  the  distance  between  the  center  0  of  the  target  and  the 

center  if  impact,  is  unbiased  with  mean  located  at  0  and  standard  deviation  of 

25  2 

a  such  that  E[e]  =  15  ft  and  VarTal  =  — ^  ft  .  If  the  weapon  impacts  at  a 
point  which  is  at  a  distance  less  than  A  ft  from  0,  the  target  is  killed  with 
probability  one.  On  the  other  hand,  if  the  weapon  impacts  at  a  point  distant 
A  *'t  from  l  (A  'he  target  is  undamaged.  The  pararei-ers  A  and  n 

estimates  and  are  equally  likely  to  be  in  the  respective  intervals  13  <  A  <  17 
ft  and  20  <  R  <  22  ft.  It  is  required  to  provide  a  two-standard  deviation 
confidence  interval  for  the  probability  of  kill  due  to  blast  PkV 

The  objective  is  to  calculate  fTp..]  and  VarTP,.!.  It  is  easy  to  verify  tha 


First  compute  U(A,B,o)  using  (19) 

U(*,fT,a)  =  j^yrxFT 

=  .417,771  (.419,243  -  .341,345) 

=  .032,543. 

Using  (20)  it  follows  that 

E[Pkb]  =l-o  U(A,R,o) 

=  1  -  (15)  (.032,543) 

=  .5119. 

To  compute  Var[Pktl],  each  of  the  variance  components  in  (26)  is  computed 
separately.  First  the  component  incorporating  Var[A]; 

i 

(-  U(A\B\a)  +  e  °  r  Var[A] 

(B-A)  (B-A) 

-  vbj? 

=  77TT15T  (,032’543)  +  ^  21-151  e  7TT  i2  $ 

=  (-  .081,358  +  .101,088 )2  (y) 

=  .000,519. 

Next, the  variance  component  involving  Var(R]  is  computed. 


13 


Pina’  I  v 


h'encp  , 


, 


(B-A) 


U(JT,B »  - 


(B-A) 


VarfR] 


-  lin) 

I-  (71-1  ST  (-032*543)  [21-15T  e  ^  M 

(.081,358  -  .062,55 2)2  (Jj-) 


.000,  118  . 


the  variance  component  involving  Var[o]  is  computed.  Thus, 


,  -  H=) 

J -  IT  e  -  fr  e 


-HI 


(B-A)  a 


}  +  U(JT,IT,7)]2  Var [o] 


1  ,  ’  0j  ‘  7^J 

rTTri^nTrey  L  5  e  ' 21  e 


4-  .  032,543  f  (^J 


2  r  2  5  > 


=  (9.097,%0  -  7.881,533)  +  .032, 543^  (ijj 


.017, 67H  . 


VarfP.h]  =  .000,519  +  .000,118  +  .017,678 


=  .018,315 


op  -  /  Var!Pkhj  =  .1353 
V  h 


pkb  'Efpkbli2«P 


kh 


.5119  +  .2706. 
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SECTION  IV 


ESTIMATION  OF  E[Pkb]  AND  VAR[Pkb]  WHEN  «x  /  a 


I .  Background 

For  general  purpose  bombs,  range  delivery  error  is  usually  greater 
than  deflection  delivery  error  a  .  From  (11),  the  expression  for  the 
probability  of  kill  due  to  blast  Pkb  is  given  by 

A  0 

n  -  —  /  re'ar  I  (hr2)  dr 


kb  a  o  -  „ 
x  y  0 


R  p 

1  t  ^B-r-,  _  -ar'  . 

+  —  f  ^  re  Io(br  5  dr 

x  y  A 


(ID 


where  If1(>)  is  the  modified  Bessel  function  of  the  first  kind  of  the  zeroth 
order. 

The  development  of  closed  form  expressions  for  E[Pkbl  and  Var[Pkb]  using 
the  result  (11)  is  hindered  by  the  appearance  of  the  integral  term  involving 
the  function  Iq(’).  However,  it  is  possible  to  obtain  an  approximate 
expression  for  Pk h  by  using  the  expanded  form  of  I^(x)  and  retaining  the  first 
three  terms.  The  resulting  expression  is  more  readily  amenable  to  numerical 
computation  than  the  original  expression  (li).  Now 


r0(x)  a  1  +  Vs?  + 


Substituting  this  last  expression  in  (11)  yields 


(27) 


1  rlir  2  -dr  , ,  hr  hr  , 

■  a  a  ta)  /  r  e  f1  +  ~~r  *  ~rr  +  •••)  dr 

x  y  A 


The  result  of  the  evaluation  of  the  integrals  in  this  expression  is  given  in 
Appendix  B.  Thus,  let 


A  2 

(A)  -  \  re“ar  dr 
1 


2  a 


(  1 


-aA 


(?«'• 


A  2 

J?(A)  =  /  e_ar  dr 

n 


_ i  r.  -a  A  /,  . /  a 

^  [  1  -  »?  ^  l  +  +  — 

a  5 


?  \ 

?  a  A  . 


( 29) 


A  Q  _ ,r2 

>)j(A)  =  f  r  e  dr 
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aA 


=  -r  [  1  -  e  i,  1  a  aA 


,  2.4  36  4.a 

,2  ,  a  A  a  A  a  A 

2  6  +  2~ 


)]  (30) 


A  2 

04(A)  =  /  r^  e'ar  dr 


“372  2a)  -  2l?'3A 


(31) 


A  ,  2 

!r(A)  =  J  r  e-ar  dr 


-  1A7  11  ./»  /~T7x  16  A  -a A 

-  „  7/7  *(A  '  ?a'  -  7  -  p 

Ha  a 


2  ’  4  n  A. 

( 1  +  "q  a  A  +  —  ■ ■  a  A  (321 


'in  '• 

/.»  r  10  -ar  , 

-( A)  =  J  r  e  dr 


16 


945  /  «  .  945  A  -aAr  f.  x  2  .2 

$  (A  /  2a)  -  "j  e  ^1  +  ^  aA 

d 


32  a 


11/12 


4  2.4  8  3.6  16  4.8-, 

+T5aA  +T05aA  +945aA)* 


(3: 


Expression  (27)  can  then  he  written  as: 


V  *  ;rrI'V4>  *pJ3<A>l 

x  y 


1  f  R 


? 

h  - 


oo  4  B 
x  y 


R-A-J  {[J1(9)  -  Jj  ( A)  ]  +  -j  (J?(ft)  -  Jz(A)j 


+  ^  ['>3^)  -  J3(A)]}  -  (^)  {[ J4(B)  -  J4(A)] 

x  y 


64 

2 


+  T  lJ5<R)  ‘  'VA)1  +  Ia  tJ6(R)  "  J6(A)U 


(34 


It  will  be  assumed  that  the  parameters  A,  B,  a  ,  and  o  are  subject  to 

x  y 


estimation  error.  Further  A,  B.  a  ,  and  o  will  be  assumed  to  be  mutual  1/ 

x  y 


independent  and  characteri zed  by  their  respective  means  A,  R,  a  ,  and  a  anc 

x  y 


their  respective  variances  Var[A],  Var[Bl,  VarTo  1,  and  VarTo  ]. 

x  y 


kecal  liny  tnat  Pkt)  in  (11)  is  a  function  of  the  four  parameters  A, 


a  ,  and  o  ,  one  may  write  (11)  as 

*  y 


Dkh(  A  ,R  ,0„  ,n„) 


x  y 


'■>  0  ' 

x  y  0 


■  re"ar  I„(br?)  dr 


,  R  n  2  . 

— ;  fa)  re  Vhr }  ,1r 

x  y  A 


'  V 


(40) 


3.  Estimation  of  Var[Pkhl 


First  expression  (36)  is  written  as 


P,.  (A,R,o  ,o  )  -  PUk(^,^,o  ,o  ) 
kb'  ’  x  ’  y  kb  ’  x  *  y 

—  3Pkb  —  3Pkb 

■  <A-»)  TT  -  *  <8-B>  TTIir-  - 

A, B,o  ,o  A, B,o  ,o 

’  x  y  x  y 

3P,h  ap 

,  —  ,  kb  .  —  ,  kb  , .  % 

+  ( a  -o  )  - -  +  ( o  -a  )  - -  .  .  '  i ) 

x  x  3 a  -r  ■=-  -  -  y  y'  3o  -r  —  — 

x  A, B,o  ,o  J  y  A, B,o  ,os 

1  x  y  J,,,x’y 

It  was  assumed  that  each  of  the  random  variables  A,  8,  o  ,  and  o  are  mutually 

x  y 

independent.  Thus,  squaring  and  taking  the  variance  on  both  sides  of  'AT 
results  in 

3P,  ?  3 P,  ,  ? 

var[pkhl  ■  (-srLff-  -  *  t-5^)  _  -  VA-tB] 

A, B, 0^,0^  A,B,ox,oy 

?  2 
3  p  3  p 

+  n-  -  Varr0x1  +  -  Var[oy]  •  :4?) 

x  A, 8,0  ,0  y  A, 8,0  ,0 

*  y  x  y 

The  details  of  the  calculation  of  the  partial  derivative  terms  appear  in 

A 1  i x  0.  The  final  •vsults  are 

A  2 

,n  a  a  P, .  -  /  re”ar  In(br  )  dr 
3 P,  .  x  y  kb  1  n  O' 

_ kb  =  _ 0 _  41\ 

3 A  ( B-A)  ox oy 

8  2 

’  re  '  K(br')  dr  -  o  o  PLh 
3P,  t  O'  x  y  kb 

*b  _  »  ... _ _ _ _ _ _ _  ul 

.'TR  (B-A)  o  o 

x  y 

j"  r?  e'ar?  (I0(br?)  -  I  (br2)]  dr 

^kb  =  .^ _ _ _ _ _  ;4h) 

3ox  2( R-A)  ol  o„ 


3 P  ‘  r'  <1r 

kh  A 


3u 


( -If.' 


2(B-A)  o  a 
x  y 


Here  !.(•)  is  the  modified  Bessel  function  of  the  first  kind  of  the  first 
order.  In  series  form  one  has 


?k  +  l 


u 


i(x)  =  y 


?k+l 


r  =  n  k ! ( k  +  1 ) !  ? 

h  ■  nay  use  the  same  technique  as  for  FfPj^]  to  evaluate  expression  (43', 
•  441 .  1  •  1,  and  ;4e'  by  using  I^(x)  and  I j ( x )  and  retaining  the  first  few 
terms,  however,  we  will  show  an  alternate  method  which  relies  on  the  Chi- 
5  : : j a r ‘jnr'ion  fur  tbp  numerical  evaluation  of  Varfp^l.  let 


’ 1  3 

n  -ar  , 

r  e  In( hr  )  dr 


v P ( ‘ ' )  rn  e-ar  I^(br^)  dr 


-hen  ( 4 3 )  through  (46)  may  be  written  as 


3P, ,  P.. 

kh  *  h 


X,(  A) 


x  y 


M,.  V11’1  M 

oB  a  o  (B-A)  "  B-A 
x  y v  1 


kh 

'  X 

kK 

y 


—r~~ - r  x  { Rd  -  X  (A)  +  Y  (»)  -  Y9(A)1 

>'  j.  ( B-A)  ' 


x  V 


- C  Y  ( R'  -  Y  (A)  -  X  !H'  *  X9;A'  1 

>  (  B-A) 


It  is  possible  to  express  (47)  and  (48)  in  terms  of 


?n 


(471 


M3) 


Mo' 


(50) 


(511 


SM 


Thr*  function  P(x^|v)  refers  to  the  chi-squared  probability  function  and  is 
tabulated  in  several  places  (see  e.g.,  [1]).  The  function  is  also  available 
as  the  PROBCHI  function  provided  by  the  SAS  computer  program  [31.  In  Appendix 
D,  it  is  shown  that 


11  2  0 
*  00  =  /  rn  e"  ar  ln(br2)  dr 
0 


r(y)r  o  (t£)‘  (n+l)(n+3)  ? 

~ n'lV  [  P(  ?alr  |n+l)  +  - r -  P( 2alr  |  n  +  5) 

2a^  ( 1 ! ) 2 


(t~)  (n+1) (n+3) (n+5) (n+7) 

_  _  _  _ 


(2!)' 


P(  2aU  1 n+9) 


(yr)  (n+l)(n+3)(n+5)(n+7)(  n+9)  ( n+1 1) 

—2 - - -  P(2alP|n+13)  +  •••  ]  (54) 

( 3! ) 2 


and 


il1))  =  /  rn  e'  '■  I.lbr7)  dr 
0  1 


(4?)  r^)  r  2 

— — [(n+1)  P(  2au  |  n+3) 
2a 


Itt)  ( n+ 1 ) ( n+3) ( n+5) 

— a - yTJV - P(  2a!)2  |  n+7) 


tM 


ITT)  ( n+1) ( n+3) ( n+5) ( n+7) ( n+g) 

— - 5T?i -  P(2aU2|n+ll)  +  •••] 


0  I  T  I 


(55) 


4.  Example 

A  general  purpose  bomb  whose  main  effect  is  blast  is  dropped  trom  air  and 
aimed  at  the  center  f)  of  a  target.  The  range  delivery  error  is  20  ft  with  a 
standard  deviation  of  3  ft.  The  deflection  delivery  error  is  17  ft  with  a 
standard  deviation  of  2  ft.  If  the  weapon  impacts  at  a  point  which  is  at  a 

distance  less  than  A  ft  from  0,  the  target  is  killed  with  probability  one. 

On  the  other  hand,  if  the  weapon  impacts  at  a  point  distance  R  ft  from  0 
'A  <  R1 ,  the  target  is  undamaged.  The  parameters  A  and  R  are  estimates  and 
are  equally  likely  to  he  in  the  respective  intervals  13  <  A  <  17  ft  and 

20  <  R  <  22  ft.  It  is  required  to  provide  a  two-standard  deviation  confidence 

interval  for  the  probability  of  kill  due  to  blast  P^. 


Let  o  and  a  be  the  respective  range  and  deflection  delivery  errors, 
x  y 


One  has 


A  =  E[A1  =  15  ft 

R  =  E[R7  =  21  ft 

j  =  Efa  1  =  20  ft  ; 

x  x  ’ 

7  =  EO  1  =  17  ft  ; 

y  y 

Calculation  of  EfP^l 

F[»kb1  is  calculated  using  (401. 


Varf  A1  ,  titm!  ,  i  ft’ 

Varm  .  mgai ,  i 

Var[oxl  =  .1  -  4  ft 

Var[oy]  =  ??  =  4  ft?. 

for  convenience,  h’l’  bar  will  on 


dropped  from  the  symbols  describing  the  input  parameters. 


From  (38)  and  (39),  a  and  b  are  computed 


*  =  j  (-7  ♦  -y)  =  j  [-4  ♦  -4l 

oj  (20  j  (17P 


=  .001,490,051,9 


4  H?  -  -?)  -  i  [-4  - 

4  ay  oj  4  (17)  ^  (20p 


=  .000,340,051 ,9 


Next,  each  of  the  quantities  J^(A)  and  J^(8),  i=l,2,...,6  is  computed  using 
expressions  (27)  through  (32) 


Computation  of  J^(A) 

2 

•I,  (A)  =  ^  (1  -  e‘  aA 


'  2)  (  .001~,497T,051 ,9) 


exp[ -( .001 ,490 ,051 ,9) ( 15)^] } 


=  95.583,609,22 


' if  a*  ion  of  ,1  ?(  A) 

*  -L  .  1  -  e‘  jA?  [i  +  aA?  +  4^)] 

a " 


.  '01 ,400,051 ,9) 


{1  -  (.715,150,923)  [1  +  .335,261,677,5 


1  1 1 


-  ’  1  -  (  .71s, 150, 923)(  1. 391  , 461  , 874)] 


4 ’'4,640, *58  x  10 


Computet i on 


•J  j  ( A) 


Conputat i on 


J4(A) 


C< input  ,t*  i .  m 

■  !S(A) 


of  J3( A) 


=  [1  -  n-  (1  ♦  aA2  +  i^+  i^l] 

<1 


12 


(.001,490,051,9)' 


[1  -  ( -715.150,Q23)( 1  +  .335,261,677, 


.11  2,400,392,4  .057,603,544 ,1  .012,633,34?,?,, 

+  *  r  +  o  A  M 


24 


12x10“ 


149,005,19) 


v-  f  l  -  (  .715,1 50, 9?3)(  1 .303,260,075)  ] 


10 


-  4.365,063,756  x  10‘ 


of  J4( A) 


(A  /2a) 


__A  -  aA2 

2a  P 


- ll - —  *r(  15)  /(2)(. 001  , 490, 051  , 9)1 

(  ?)  (  .091  ,490 ,051  ,9) ^ 

f  1 5 ) ( .'’15,150,923) 

rn  'i  '  rvM  r.  in  OS  t  4~ 

V  •  *  *  1  '  *  »  7  ■ 


(/ti  )  ( .293,57) 

( 2) ( .000,057,517 ,3) 

‘•23.661  ,527,2 


3,599.627,505 


1  5/  r 


7/2 


:(A  /2a) 


15 

B 


a2  AA ) 


24 


(8) ( .001  ,490,051 ,9) 


jj2  *( .818,854,904,7) 


-  i  [1  +  ^  (.335,251,677,5) 

(8) ( .001 ,490 ,051 ,9) J  J 

+  ( .335,261 ,677,5) 21 

=  7.639,817,718  x  10Q  -  ( 6.079,754 ,57) ( 1 .253 ,481 ,223)  x  109 


=  1.895,952,4  x  10' 


Computation  of  ,J,(A) 
n 


,  tA\  _  945  /it  wa  /—  945  A  aA  r .  .  ?  _*? 

J6(A)  "  ~  \\j?.  ^(A  /2a)  ‘l2Te  [1  +  3  aA 

oc  a  a 


4  2.4  8  3.6  16  4.8, 

+  !5aA  +T05aA  +  945  3  A  ) 


_ 945  /n 

(32) ( .001 ,490,051 ,9) 


Y3/2  ♦(. 818, 854, 904, 7) 


( 32) ( .001 ,490,051 ,9) 


Tp-  [1  +  |  (  .335,261  ,677  ,5) 


+  yy  (.  11  2,400,392,4)  +  y— -  (.037  ,683,544,1  ) 


+  qjl  (.012,633,848,2)] 


(.293,57) 

( 32) ( .001  ,490,051 ,9) 1 ^ 

(32 )( .001, 490, 051, 9)b 


J,(A)  =  54.195,192,55  x  10lS  -  54.193,750,53  x  101S 

o 

=  1.442,016,6  x  1012  . 


(21) (.51 8,. 145, 698, 4) 
(2)( .001  ,490,051 ,9) 


J4(b) 


=  (T)~(~.000,057 ,517 ,8)  -  3,652.644,47 

=  5,765.317,71  -  3,652.644,47 
=  2,112.673,24  . 


Computation  of  JC(B) 

_ 3 

i  la\  13  /tT  a,d  .  15  R  -  aB?  2  „2  4  2  4. 

■1^(8)  =  — y/2  $(B^2a)--^-^e  (l+--aR  *  ^  al  ) 

8d  d 


15  'I  TT 


8(.001,490,051,9)7/2 


$(1,146,196,867) 


(.001,490,051,9) 


[1  +  {  (  .657,112  ,887,91 


+  tt  (.431,797,147,4)] 


=  9.737,599,188  x  109  -  9.582,287,423  x  109 


15.531,176,54  x  10  . 


Computation  of  J  (B) 
o 


J6(B)  =  ~~~ll72  ${R  /2a) 
32  a11' * 


945  B  -  aR  ; .2  d2 
~ I?  ~5  6  [  1  +  aR 

a 


+  _i  aV  +  _«  16  JL6  4  8- 
+  15  3  R  +  105  3  R  *  945  3  R  J 


_ 945  /n _ 

(32)(. 001, 490,061,9)' 


T  $(1,146,196,8671 


(.001,490,051,9)' 


( .618,145,698,4) 


1  +  4  (  .657,112,887  ,9)  +  t£  ( .431  ,797 ,347  ,4) 


+  ln^  (.283.739.60?)  +  q-^  (.186,448.949.3)] 


=  69.076,394  .55  x  1015  -  (43.763  .661  .05)  (  101S)  (  1.577.996.7.81) 


-  (69.076.394.55  -  69.058.894.38)  x  10 


=  17.500,174.5  x  10l  . 


co 'i!  ite  rr P,.s] ,  the  numerical  values  of  ;  .  j  .  h,  ,)  (4)  an 6 
■  ■’  x  y  1 

J^(B)  (i  =  l . 6)  are  substituted  in  (40).  To  simplify  the  computation, 

the  three  major  terms  of  (40)  are  computed  separately. 


T  e  rm  1  = 


+“4  J?(A)  +'6“4  J3(A)] 


_ J _ [or  ran  rnq  O?  +  'l- — i?3Sl  .051  ,  9 ) 

C 17)  ( 20)  +  4 


(  1.479,540,658)  x  106  +  (4.365,863,756)  x  10 


r 95.583.609.??  +  .n?\,3i4.60i  +  .non,nn?,?65,?l 


=  .281,19  . 


T-rm  ?  =  (rRx)  {[  J  (B)  -  Jj(4)  ]  +  ^  [J?(B)  -  J?(A)] 

'  <  y 


*  ^  i J 3 ( R )  -  J3(A)j: 


n7)}?0)  (  ?T- 15 )  [  ( 161.623 ,3-30.6  -  95.583,609,??) 


+  »•-" x-0-5-1-^.  (8. 80S, 900, 503  -  1  .478,540.658)  x  106 


+  — 4^°51,q^  (96. 896. 017. 28  -  4.365,863.756)  x  1010] 


rrrrfi^Tfii 


>.039,721 ,38  ♦ 


+  ( 2. 400, 5 19) 4 (92. 530, 153, 52) j 
(64)(106) 


(21) 

(17)  (20) (6) 


(66.039,721,38  +  .105,545,315,3  *  .000.048,00 


( 21 ) ( 66 . 145 ,314,59) 


.68091 


Term  3  =  (g^)  {[ J4(B)  -  J4(A)]  +  ~  [ J5(B) 

x  y 


Jc(M] 


&  K(B)  ' 


flTTITO)'  UT7ib)  [  (2 ,112.673,24  -  923.661  ,527  ,2) 


+  (  15.531  ,176,54  -  1.895,952,4)  x  107 


+  (  17.500,174,5  -  1.442 ,016 ,61  x  1012] 


_ 1 _ r  i  1Rq  011  7n  +(?^,519),ln^l8.224,H) 

(  17)  (  20)  (6)  t  1,189.011.713  (  4 )  ( 10) 


+  (_2.400 L5 V) )  (  16.058,157,9)  j 
(64)(104) 


7TVw4nwV\-  (1.189.011,713  *  1.964.3?!  ,57  +  .000, 833, 175. ^) 

(17)1  1(D) 


1.190.976.868  _  o, 

U7T(ToT(6)“-  -Sm3'81  * 


E.1  P.  .  ]  =  Term  1  +  Term  2  -  Term  3 
kr> 

=  .281,19  +  .680,91  -  .883,81 
=  .378.29  . 

0a  1  cu  1  a'  ion  of  Varr 1 

kh 

To  compute  Varf  Pkhl,  the  following  procedure  is  used: 
a.  Retaining  only  the  first  two  terms  of  the  series  (84)  and  (85)  one 
finds  f'-on  the  tables  (e.g..  Til)  using  interpolation  the  following 

P(2aA2|21  =  0.284,843,021,8 
P(2aA213)  =  0.119,900,605,6 
P(2aA?|5)  =  0.015,452,739.7 
P(2aA2|6)  =  0.004,893,986,9 
P(2aA?|7)  =  0.001,442,416,9 
pf?o.a?  o'-  =  fi.nnn,i09,195,8 

Similarly,  one  finds 

P(2aR2|2)  =  0.481,647,535,2 
p( 2aR2 1  3)  =  0.274,248,026,6 
P(2aR2|5)  =  0.066,542,951,0 
P(?aR2l6)  =  0.029, 130, 319, 2 
P( 2aR2 | 7 )  =  0.011,944,485,8 
P( 2aB2 | 9)  =  0.000,109,195.8 


M 


b.  Substituting  these  values  in  (54)  and  (55)  one  obtains  for  n=l,2  the 


f ol 1 owi ng : 


Xj ( A)  =  95.602,888,58 
Xj ( B)  =  161.747,910,3 


X?(A)  =  923.977,344,6 


X?(B)  =  2,115.029,112 


Y?(A)  =  14.428,854,99 


Y?(B)  =  61.985,977,35 


c.  These  values  are  next  substituted  in  (49),  (50),  (51),  and  (52' 


yi el d : 


£-  =  0.016,176,955 


=  0.016,248,075,3 


-0.015,179,030,5 


-  =  -0.016,486,370,3 


4.  Using  now  (42)  results  in 


VarT|>  1  --  (0.016, 175  , 955)?  (4)  +  ( 0.016 ,248 ,078 , 5)  ?  (~) 


V  (-0.015, 170, 030, 6)?  (9)  +  (-0.016 ,486 ,370,3) 


0.003.597,710,3 


...  t  *.  «*.  V.  -*\ 


SECTION  V 


CONCLUSION 

In  this  report,  the  probability  of  kill  due  to  blast  of  an  exploding 
weapon  is  related  to  the  aiming  error  for  the  weapon  in  order  to  obtain  a 
closed  form  expression  for  the  net  probability  of  kill  due  to  blast,  Pkh. 

Using  the  Taylor's  series  estimation  procedure,  expressions  for  ErPkh] 
and  VarrPk|)]  are  obtained  in  order  to  arrive  at  two  standard  deviation 
confidence  levels  on  Since  closed  form  expressions  for  ETp^]  and 

Var[Pkbl  are  not  possible  to  derive,  several  approximation  techniques  are  r 
to  obtain  numerical  results. 
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APPENDIX  A 


EVALUATION  OF  K( r) 


The  intent  is  to  evaluate  the  integral 


2n 


Hr).  I  .  &iA)]  d6 


0 


?a 

x 


?.a 


Using  the  trigonometric  formulas 


?  1 

COS  '9  =  J  (1  +  cos  29 ) 


sin29  =  j  (1  -  cos  29) 


gives  for  K( r) 


Kfr)  =  exp[~  (-£  +  -ij]  /'  exp[~  (-J-  -  -~)  cos  ?o]  do 


a"  a"  0 

x  y 


a  a 

y  x 


,  .  1,1  1-1 

Let  a  =  ^  (—  +  — ) 

o'  o' 
x  y 


b  =  4  (T  -  ~2}  • 


a  a 

y  x 


Then 


K(r)  -  e 


2  ?n  u  ? 
ar  t  hr  ‘cos  29 


dO 


in  K(r),  let  2o  -  x,  then 


APPENDIX  B 


EVALUATION  OF  .^(A)  THROUGH  Jfi(A) 


In  this  appendix,  expression  (26)  which  is  rewritten  is  to  be  evaluated. 


1  ,A  ar2  , ,  .  b2r4  .  bV 


kb  “  a  o  . 
x  y  0 


r  -  ai  /  i  u  r  u  r  ,  \  , 

J  re  (It  — j-  +  +  •••)  dr 


1  ,  1  .  “  -  ar2  „  b2r4  b4rS  ,  „ 

+  ~  ^r-t^  >  re  ^  +  ~  +  ^r  +  dr 

x  y  A 


.  ,  B  „  2  ,24  48 

“  o  0  J  r  e  ^  +  4“  +  “64“  +  '•*)  dr  •  (R-0 

x  y  A 

It  is  clear  that  this  expression  is  reducible  to  the  evaluation  of  six 
integrals,  namely 

A  2  A  2 

J, ( A)  =  /  re"  r  dr  ;  J,(A)  =  /  rb  e“  ar  dr 

1  0  ^  0 


A  9  -ar2 


A  2  -ar2 


Jo(A)  =  /  re"  dr  ;  J.(A)  =  /  re”  dr 


i  r  6  -  ar2  ,  1  ...  r  10  -  ar2  , 

Jc(A)  =  /  r  e  dr  ;  ( A )  =  J  r  e  dr  . 


0 


0 


Each  of  these  integrals  will  be  evaluated  separately. 


A  _  2 

1.  Evaluation  of  J,(A)  =  /  re”  ar  dr 

1  0 


J  ( A)  =  /  re 
1  0 


A  2 

ar 


dr  =  -  J-  JA  die" 
2a  n 


1  -  ar' 

2a  ’’ 


l  .  -aA' 

21  ^  -  e 


B-  , 


37 


40 


J,(*>  =  +  2T  C-  2^)  /  d(e  -  ) 


,5  „2  2  ,  A  A  2 

*r  e-  aA  -  -5=  [  e-  ar  r3  -  /  3  r2  e"  ar  dr 

4a'  0  0 

a!  -  aA2  5_  -  aA2  15  rA  2  -  ar2 

2ae  -4  j2«  *4a2/0' 

5  2  12  A 

A  -  aA  5  A  -  aA  15  ,  K  ,  -  ar 


A^  -  a A2  5  A^  -  aA? 

2a  0  "4  2  e 

a 


,c  2  A  A  2 

3  L  e  r  -  /  e  dr] 
8a  0  0 


5  2  2  2  2  A 

An  -  aA'  5AJ  -  aA'  15A  -  aA'  x  15  ,  -  a 

~7a  6  - 2  6  "  ~  e  +~ Je 

'  4a'  8aJ  8a1  0 


Making  the  change  in  variable  /a  r  =  u//2  in  the  last  integral  yields 


,  A5  -  aA2  5A3  -aA2  15A  -  aA2  x  15/ n  1  / 

J5  A  =  ~  Ta  e  -7Te  ~7“3e 

4a  8a  8a  /2n  0 


ising  e<prossion  (8-5)  one  obtains 


1  (h\  -  15/ rr  15  _A  -  aA  r.  2  ,2  __4  2.4-, 

J5  A  "  q  7/2  ^A/2a^  "83  (1  +  3  aA  +  15  a  A  ]. 

o3  3 


A  2 

>.  Evaluation  of  J  (A)  =  j  r^3  e"  ar  dr 


)fi(M  =  J  pl;)  e_  ar  (1r  *  -  (,yj)  /  r<1  <ke"  dr  ) 


APPENDIX  C 


EVALUATION  OF  THE  DERIVATIVES  OF  Pkb 

The  purpose  of  this  appendix  is  to  derive  expressions  (49)  to  (5?). 
For  convenience,  these  expressions  are  repeated  below.  Let 


U  n  2  , 

x  ( U)  =  /  rn  e"  ar  In  (br2)  dr 

n 

u  ?.  , 

Yn(U)  =  /  rn  e”  ar  I.  ( br  )  dr  . 

'  r\  ^ 


(C-l) 


(C-  ?) 


Then,  expressions  (49)  to  (52)  which  are  to  be  proved  are 


3  A  B-A  a  a  (B-A) 

x  y 

!V  *i(B)  pki> 

an  '  a  a  (R-A)  " 

x  y 


(49) 


(50) 


3  P 


kb 


1 


3a 


o  l  "o 

2a  a  (B-A)  * 

x  y 


[X?(B)  -  X- ( A)  +  Y  (B)  -  Y  (A)] 


3  P 


kb 


1 


[ Y-( B)  -  Y?( A)  -  X„(B)  +  X„( A) ] 


9  y  2axaZ  (B-A)  2222 


(51) 


(52) 


The  derivation  is  divided  into  five  sections.  In  Section  I,  some  recurrence 
relations  are  derived  which  will  he  found  useful  in  Sections  IV  and  V. 
Section  II  to  V  provide,  respectively,  the  derivations  of  (49),  (50),  (51), 
and  ( 52) . 


Section  1.  Some  recurrence  relations 

Using  Leibniz's  rule  and  applying  it  to  (C-l)  it  is  found  that 


| 


-  -  /  r  e  I  (hr') 

n 


The  integral  portions  of  (C-3)  is  recognized  as  being  xn+?  (D)  ,  hence 


ax  ( u) 

n 

3  d 


xn+?(l))  • 


In  a  similar  fashion  one  may  differentiate  (C-l)  with  respect  to  b  to  obfain 


ax  (M)  II  2 

" _  n+2  -  ar  , 


j  r"Tt  e"  I . ( hr  )  dr  . 


This  co’vs  about  owing  to  the  manner  in  which  Iq($(x)1  is  differentiated, 
namely , 


dInU(x)) 


\  V  \  n  I  I  J  .  /  „  \ 

dx  dx  * 


The  final  result  is 


a/  (d) 

-A  ...  =  Y  (U) 
ab  n+2  U 


(C-7) 


I h (»  definitions  of  a  and  b  that  are  given  by  (Q)  and  (10)  may  readily  be 

rearranged  to  u’V’ 

2( a-b)  =  ~  (C-R) 


?( a+b)  =  -j  . 

a 

y 


Appropriate  differentiation  of  (Q)  and  (in)  yields 


Some  essential  identities  are  now  generated.  The  rule  for  integration  by 


parts  is  given  by 


/  u  dv  =  uv  -  /  v  du 


(C-14) 


In  (Cl)  let 


and 


rndr  =  du 

-  ar2  ,  2, 

e  1q( br  )  =  u 


( C— 1 5 ) 


(C-1M 


It  foil ows  that 


n+1 


n+1 


and 


2  2 

du  =  [ 2br  e"  ar  I^br2)  -  2ar  e"  ar  IQ( hr2) ]  dr  . 


tC-17) 


(C-H) 


Substitution  into  (C14)  gives,  upon  imposing  the  limits  given  in  (C-l)  : 

,,n+l  -  aU2  In{bll2)  _ 

<„<">  ■ - - *  nSf  W'  -  nM  ’n*?"1'  '  <C-’q> 


A  multiplication  by  (n+1)  then  gives 


(n+1)  Xn(U)  =  Un+1  e'  aU  IQ(  bll^)  +  ?aXn^2(U)  -  ?bYn+?(l!) .  (C-'>'» 


By  taking 


r  dr  =  dv 


(C-2 1) 


e"  ar  Ij(br2)  =  u 


( C-22) 


one  may  also  integrate  (C-2)  by  parts.  In  this  case 


( C-23) 


?  2 

du  =  [  2br  e”  dr  I^br2)  -  e~  ar  Ij(br2)(2-  +  2ar)]  dr. 


(C-24) 


The  latter  relation  follows  since 


dIiU(x))  IA*(*))  . 


(C-2 5) 


After  minor  algebraic  manipulation  the  integration  by  parts  gives 


( n-1)  Yn(U)  =  Un+1  e'  aU  I^bll2)  +  2aYn+2(U)  -  2bXn+?(U)  .  (C-2fi) 


The  sum  of  (C-20)  and  (C-26)  yields 


Z(a-b)  [Yn+2(U)  +  Xn+2(U) ]  =  (n+1)  XR(U) 

2 

+  (n-1)  Yn(U)  -  Un+1  e*  aU  [  IQ( bU2)  +  I^blJ2)] 


(C-27) 


while  their  difference  yields 


?(^b)  [Yn+?(U)  -  Xn+2(U)]  =  (n-1)  Yn(.n 

? 

-  (n+1)  Xn(U)  +  Hn+1  e*  aU  [ln(bU?)  -  I ( bnpT ]  .  (C2«) 

Identities  (C-27)  and  (C-28)  will  prove  to  be  important  in  computing  (51)  and  (521 


46 


■  '■  y  r. 


Section  II. 


3Pkh 

Computation  of  — 


By  incorporating  ( C- 1 )  into  (11)  one  obtains 


BXj(B)  -  AXj(A)  -  X ^ ( B)  +  X?(A) 


kb 


(B-A)  o  a 
x  y 


(C-29) 


EtPkbl  m  Section  I V- 2  may  also  be  calculated  using  (C-29).  Rearranging  one  obtains 


(B-A)  P 


BX1(B)  -  AXt(A)  -  X2(B)  +  X?(A) 


kb 


a  a 
x  y 


(C-30) 


Now  differentiate  (C-30)  with  respect  to  A.  Proceeding  formally  one  gets 


(B-A) 


3P 


kb 


3 A  ‘  Pkb 


ax  (A)  ax  (A) 

-A  -  X  j  ( A)  + 


3  A 


3  A 


a  a 

x  y 


( C-31 ) 


From  Leibniz'  rule  one  gets 


and 


3X1(A)  -aA"  2 

3A  =  A  *  J0(bA  > 


3X?(A)  2  -  aA2  2 

3A  =  A2  e  !0(bA2) 


( C- 32 ) 


( C-33) 


Substitution  into  (C-31)  gives 


3  P , 


kb  *i(s) 

(»-*>  7T  '  >“kb  *  '  sV 

x  y 


( C-34) 


or  in  rearranged  form 


3  P ,  ,  j  j  P .  -  X .  ( A ) 

_ kb  _  _x__y  _kb  __  1 

3A  (B-A)  o  o 

v  x  y 


(C-35) 


which  is  the  same  as  (49) 
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-  'X 


Section  III. 


Computation  of  — r-„— 

o  l) 


Returning  to  (C-30)  and  differentiating  with  respect  to  R  gives 


3P,  , 

^B'A^  HT  +  Pkb 


3X.(B) 

8  +  xi(B) 


3X2(R) 


(C-36) 


Since  by  Leibniz'  rule 


jX,(B)  ? 

,  Be’  a  In(bR2) 


(C-37) 


3X?(R)  R2e‘  -«2  T  ,KR?, 
— gjj~  =  R  P  I n(  bB  ) 


(C-3R) 


one  finds  that 


3P.  K  X, (B)  -  a  j  P, . 
kh  1'  x  _y  kh 

"3R'  '  "  "( B-A")  o  a  ' 
x  y 


( C- 39 ) 


This  is  the  same  as  (50) . 


Section  iv. 


Computation  of 

x 


The  computation  of  - — -  is  more  involved.  Start  with  (C-30)  and 

o  0 

X 

differentiate  to  get 


(B-A)  r-k-h-  =  rRX.(R)  -  AX.  (A)  -  X„(B)  +  X_(A)1 

'10  O  0  do  1  I  /  ? 

x  x  y  x 


-  -y--  [BXj  ( B)  -  AXt  (A)  -  X2(R)  +  X?(A)] 


(C-AO) 


By  the  chain  rule 
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!i‘_!  _  3J •_!  It)  9b 

8ox  3a  3ox  +  3b  3a^ 


and  owing  to  (C-4),  (C-7),  (C-10),  and  { C-l 1 )  one  gets 


»xn(u)  . 

-55- "73  iV2<u>  *  V2(U» 

X  ,/0x 


A  similar  calculation  gives 


3Xn(U)  1 

~Ta  =  n  fXn+2^IJ*  ‘  Yn+?(U^  * 


y  2a' 


The  complicated  substitution  into  (C-40)  gives 


(B-A)  — ^  =  — J—  [B(X3(B)  +  Y3(B)) 

*  y 


-  A( X„ ( A)  +  Y,( A) )  -  (X.(B)  +  Y.(B))  +  (X,(A)  +  Y,(A);] 


(B-A)  PL 


Multiplying  by  -g— ^  and  then  invoking  (C-8)  gives 


d  P 


[B(X3(B)  +  Y3(B) )  -  A(X3(A)  +  Y3(A))  -  (X4(B)  +  Y4(B)j 


♦  CX.(A)  +  Y,(A))1 


( C-4] ) 


(C- 42) 


rc-43) 


( C-44 ) 


(C-45) 


With  n  set  to  1  and  then  2  in  (C-27)  and  these  results  inserted  into  (C-45) 
one  gets  after  some  nice  cancellation 
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Replacing  Pkh  by  (C-29)  yields  after  cancellation 


3P 


kb 


x  3a  "  2( B-A)  a  a  1  A 2 
x  x  y 


[-  X  (8)  -  Y  (B)  +  X  (A)  +  Y  (A)]  .  (C- 


Di vision  by  then  gives 


3Pkb  [ X2( B)  -  X2<A>]  +  t Y?{B)  -  Y2(A)] 


3a 


2( B-A)  a‘ay 


(C- 


This  is  the  same  as  (51) . 


3P 


Section  V. 


Computation  of 


kb 


3o 


3P 


kb 

The  computation  of  ^ —  follows  along  the  same  lines,  but  uses  ( C-43) 


(C-9),  and  (C-28)  rather  than  (C-42)  ,  (C-8),  and  (C-27).  The  result  is 


aPkb  ^X2(P)  ‘  X2(A)]  ‘  ^Y2(B)  '  Y2(A)^ 


3a 


2 ( B-A )  o%o‘ 


(C-< 


which  is  the  same  as  (52). 


APPENDIX  D 


DEVELOPMENT  OF  EXPRESSIONS  FOR  Xm(U)  AND  Ym(U). 

In  this  appendix,  relations  (54)  and  (55)  are  derived.  First  expression 
(54)  is  developed.  Define 


U  2  0 

Xp(lJ)  =  /  rf1  e"  ar  IQ  (brZ)  dr 


(0-1) 


and  make  the  change  in  variable 


p  =  ar  ; 


(D-2) 


this  gives 


aU 


2  n-1 


VU>  =  - P  2  e-p  Iq(J  p)  dp 

2a  2 


( 0-3) 


The  Maclaurin  series  expansion  for  Ig  (—  p)  is 


,  -  .  2k  2k 

>„c!  »>  ■  r  (t)  — ' 1 


k=0  a'  ( k! ) 222k  * 


(0-4) 


Substitution  into  (D-3)  followed  by  interchange  of  the  order  of  integration 
and  summation  gives 


xn(u) 


1 


l 


(!) 


2k 


air 


-  p 


4k+n-l 

- 


2a 


— ^  k=0  (k!)Z2Zk 


dp 


(0-5) 


The  integral  portion  of  (0-5)  is  recognized  as  an  incomplete  gamma  function 
and  may  be  evaluated  using  the  tables  of  the  Chi-Squared  Probability  function 
expressed  in  (53).  To  see  this,  consider  the  integral 
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(0-6) 


p(a,x)  =  /  e*  P  p01"1  dp  * 

Let 

x  =  (0-7) 

and  a  =  ^  .  (0-8) 

This  allows  (D-6)  to  be  written  as 

2 

X_  v  i 

P((£.  4)  =  —“/  ?  e'  P  P2  ’  dp  •  (°“9) 

r(f)  0 

Upon  replacing  p  by  ^  yields 

2  v2  -  i  f  -  1 

P(f.  ~ - /  e  2  W  dt 

?  0 

2  r(|) 

=  P(X2|v)  (0-10) 


as  defined  in  (53).  Hence  (0-6)  may  be  written  as 
P(a,x)  =  P(?x|2a)  . 

Comparison  of  (0-6)  with  (0-5)  indicates  that 


2a  = 

4k+n+l 

(0-11) 

and 

2x  = 

2 

2alJ  . 

(0-12) 

Hence  (0-5)  may  be  written  as 
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-  r(^*in±I)  P(2aU2| 4k+n+l) 

k=0  (k!)2  22k 

But  owing  to  the  properties  of  the  gamma  function 

r4k+n+K  _  (n+l)(n+3)  •••  (n+4k-l)  ,n+l, 

ri  2  >  =  .,2k  ri  2  J 

for  k=l ,2 , . . . ;  and,  hence. 


xn(u)  - 


1 


2a 


n+1 

2 


(n-13) 


(D-14) 


Xn(U)  =  “n+T  [p(2alj;?ln+1) 
2a~ 


^  I 

k=l 


l4aJ 


2k 


(n+l)(n+3)  •••  (4k+n-l) 


(k!)' 


P(2air|4k+n+l)  ]  . 


(H-15) 


This  is  the  same  as  (54) . 

The  development  of  (55)  follows  along  the  same  lines  except  that  in  (D-3) 
Xn(U)  is  replaced  by  Yn(u),  as  defined  in  (48),  and  IQ(j  p)  is  replaced  by 
p)  where 


hCf  p)  -  l 

1  k=0 


2k+l 

_ e _ . 

(k) !  (k+l) !  22k+1 


(0-15) 


The  final  result  is  the  same  as  (55)  or 

.  2k+l 

2  '  ;  (.;) 


V„(U)  . 


r(^) 


n+1 


l 


(n+1 ) ( n+3) • •• (4k+n+l ) 


k=n 


k!  (k+l)! 


P(2air|4k+n+3)  ( T>-17) 


2a 
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(The  reverse  of  this  page  is  blank.) 


